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A vortex-glass transition due to point disorder in layered superconductors is studied for the 
case with an applied field parallel to the layers. Our calculation of tilt responses indicates that, 
irrespective of the magnitude of the field, the resulting glass phase, Josephson-vortex-glass 
(JG), should have a transverse Meissner effect, as in a planar splayed glass phase, only for a 
tilt perpendicular to the layers. Further, focusing on the high field (and/or high anisotropy) 
region B\/T > (po/d?, where F is the mass anisotropy in the Lawrence-Doniach model, the 
JG transition line Tjq{B) is shown to have a similar form to a B-T line following from the 
disorder-free Lindemann criterion and to decrease with increasing By/V, in marked contrast to 
the disorder- free melting line insensitive to bVt in such the high field region. This Tjq{B) line 
seems to have been recently observed in a.c. susceptibility and in-plane resistivity measurements 
in BSCCO and qualitatively explains a field dependence at lower temperatures of previous 
BSCCO resistivity data showing the so-called in-plane Lorentz force-free behavior. 
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§1. Introduction 

In bulk type II superconductors, different types of superconducting glass phases with no Ohmic 
dissipation are expected, depending on types and configurations of the static disorder (randomness), 
to be realized in a magnetic field. At least in homogeneous materials (i.e., with no granularity), 
no glass ordering is expected in zero field without transition to the Meissner phase because there 
is no origin of frustration. By contrast, in a nonzero field a spatial variation of |^| due to field- 
induced vortices makes a formation of a glass phase with vanishing resistivity rather favorable at 
a low enough temperature.BH& In this sense, the superconducting glass phases are peculiar to the 
nonzero field case and hence, can be called the vortex-glass phases.0'& 

So far, Ginzburg-Landau (GL) studied' of the vortex-glass transitions in a layered supercon- 
ductor were limited to the most familiar case (B || c) with a field perpendicular to the layers. 
No GL description of the case (B _L c) in a field parallel to the layers and with no disorder was 
available even semiquantitatively until recently, and it was first performed in ref.5 in an approxi- 
mated but systematic manner focusing on the lowest Landau level (LLL) modes of the pair-field ■i/' 



(superconducting order parameter). Since a systematic description of the glass transitions in B || c 
is available at present based only on the LLL approach, it is natural to, based on the results in 
ref.5, extend our study of a vortex-glass transition to B _L c case. It is our purpose in this paper 
to perform this. 

The present study was motivated in part by resistive data of high Tc materials. First, the 
vanishing behaviors of Ohmic dissipation in real YBCO materials in B _L c are controversial. 
Resistive data& of 90 K YBCO for a current perpendicular to B and parallel to the layers have 
indicated a continuous thermodynamic transition, while the c-axis resistive (pc) data in ref.7 have 
shown a first order transition as well as in B || c. More recently, a resistive measurement^^ in 
YBa2Cu408 has shown an enhancement of a sharp vanishing, suggestive of a thermodynamic first 
order transition, of pc with decreasing the tilt angle between the B direction and the CuO planes. 
Such a vanishing, at least, of pc in B _L c cannot be understood without point disorder on the CuO 
layers.^ Further, since the vanishing of resistivities in vortex systems is generally of a collective 
nature, a collective disorder effect near a freezing transition must appear irrespective of the current 
direction. Clearly, these observations!' Bi^ need to be discussed in the context of a point-disorder- 
induced glass transition. 

Second, to our knowledge, the issue of the so-called in-plane Lorentz force free behavior0) in 
BSCCO in tesla range has not been resolved sufficiently. In tesla range of this material satisfying 
see §2), the resistive data are quantitatively independent of the relative angle 
between B and the applied current, both of which are applied in directions perpendicular to the c 
axis. This "force free" behavior was verified later in refs.ll and 12. To explain the data in ref.lO, 
a (technically inevitable) deviation of the applied field from the CuO plane© and the presence 
of the so-called double-kink excitation© were argued in literatures. The formeiS-' interpretation 
ascribing the angular independence to a vortex flow of deviation-induced pancake vortices inevitably 
contradicts the tendency seen in Fig. 2 of ref.ll that, at higher temperatures (around 86K) below 
Tco, the field dependence of resistivity disappears with increasing B. Further, the picture, used 
in ref.l2, of a vortex liquid confined within the interlayer spacings with double-kinks thermally- 
excited canno explain Ohmic dissipation seenS'© in the disordered phase. In ref.l5, one 
of us has argued based on the anisotropic GL fluctuation theory that, with no static disorder, 
the in-plane Lorentz force free behavior is due merely to a strong enough anisotropy and will 
simultaneously be accompanied by the absence of a i3 dependence of resistive data in the same B 
range. Although a statement similar to this remains valid in the Lawrence-Doniach (LD) model 
with layer structure inducing the intrinsic pinning effect (see §2) and explains the above-mentioned 
behavior (B-independent "force free" behavior) in ref.ll at higher temperatures, the resistive data 
in Fig. 2 of refs.ll and 12 seem to show a B-dependence at lower temperatures with the "force 
free" phenomenon kept. Within the LD model, this i?-dependence will be clarified only by taking 
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account of point disorder on the CuO planes inevitably contributing at lower temperatures. 

A vortex-glass transition temperature is determined as the temperature at which the uniform 
glass susceptibility XG(k = 0) diverges on cooling, where xgO^) is the Fourier transformation of 
the correlation functionB^ 



Gcimd, R±) = dJ2 J d^/dzj < ^p*{r±) i)j+^{r_L + R^) > 



(1.1) 



defined consistently with the Lawrence-Doniach model 



Wld = ^E y"dyd^ 
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(1.2) 
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where d is the spacing between the superconducting layers lying in y-z plane, ipj{T±) with r_|_ = (y, z) 
in the first representation of eq.(1.2) is the pair-field (superconducting order parameter) on the j-th 
superconducting layer, eq — (T — Tco)/Tco, T the mass anisotropy, V'(r) a spatially continuous pair- 
field in the second representation of eq.(1.2), and b^^'^ is written, in terms of the zero temperature 
penetration depth A(0) or of the 2 D Ginzburg number (i.e., fluctuation strength) , in the form 
32tt^ X"^ (0) / ipQ = 2t: de^Q /k-Q,Tco. Further, in the first representation we put A-,- = Bjdy + 6A{jd,r±) 
by assuming B || I _L c (= x), where 6 A is an external gauge disturbance necessary in defining 
linear response quantities. In obtaining the second representation, the Poisson summation formula 
was merely used with no approximation. Hence, the Josephson vortices with no singular cores 
in the first representation are expressed, for convenience, in the second representation as singular 
vortices of il^{r). In ref.5, a field-tuned sequence of structural (mean field) transitions between 
various disorder-free vortex solids was examined in terms of the second representation. 

Point disorder on the superconducting planes is introduced by adding, e.g., a random scalar 
potential term0) to eq.(1.2) 



firp = t^E y"dydznj(r_L) |V'j(rj 



(1.3) 



where the random potential uj has zero mean and Gaussian ensemble Uj{r±i) ui{r±2) = 
S^qA 5^'^\r±i — r±2)5j^i. Here we did not introduce other uncorrelated disorder terms because, as 
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will be seen later, our main results in this paper are independent of the presence of such additional 
terms. 

In §2, the model (1.2) with the term (1.3) is examined consistently with ref.5 using the LLL 
approximation, and the glass susceptibility XG'(k) is examined to find a glass transition line of the 
Josephson vortex system, i.e., of B _L c case of the LD model by focusing on the nontrivial high field 
regime B > 0o/(27r(i^-v/r)- In §3, it is demonstrated by examining tilt responses of the LD model 
and of the closely related uniaxially periodic modeill' that the resulting Josephson-vortex-glass 
(JG) phase should have a transverse Meissner effect (TME) for a tilt perpendicular to the layers, 
while the response for a tilt within the layers (perpendicular to c axis) remains nonsuperconducting. 
As intuitively expected by identifying the superconducting layer in the present case with a well- 
defined splayed plane^^ formed by line-like defects, such an anisotropy of tilt responses is similar 
to that in the planar (Gaussian) splayed-glass phase examined in ref.4, although they are different 
in the origins of TME from each other. In §4, the resulting B _L c phase diagrams to be realized 
in a manner depending on the material parameters are discussed in details, and consistencies with 
available experimental data0'lli^ are pointed out. 

§2. Glass Transition Line 

Let us begin with rewritingP^ the LD model (1.2) (with no magnetic screening) by assuming the 
continuous pair-field ip{r) to be composed of the LLL modes. The disorder terms, eq.(1.3), will 
be omitted for the moment. We will express tp as ^(r) = J2Q,q^ fi^j ^±)'^Q{xiy) ^wi^Qzz) with 
the LLL eigenfunction UQ{x,y) oc exp{iQy — \/r(x + (5r^)/(2r^)) in the gauge A = Bxy and 
substitute it into the second representation of eq.(1.2). Due to the layer structure with period 
d, the momentum Q is writteni) as Q = r'^dn + qy with = ^0o/27ri?, where \qy\ < dr^^/2, 
and the partial LLL degeneracy is measured by the integer n with the degree of degeneracy = 
Lx/d, where denotes the linear system size in /i-direction. For the moment, a constant gauge 
disturbance 6A will be assumed. Using the Poisson summation formula, the quadratic part of 
eq.(1.2) is expressed as 

+ ^ ^1 - e-P/^ < exp(-p(m + 1/2 - qyr^/df) + f -qz\ ]\ vin, qx) |^ (2.1) 

where 

Vfd^ 27rd^VfB 
p=^— = (2.2 

is a dimensioless field playing important roles in B _L c, 93(72, q^) was rescaled in the way 
Lp/^h{qy) LP, h'^ily) = Emexp(-p(m - qyT^/dY), and < s{m) >m denotes h{qy) x 
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X^mi^l"^) exp(— p(m — qyT"^/ d)"^)]. Assuming e <C 1 and neglecting higher order terms in 
qy, eq.(2.1) is simphfied in the form 
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Q 




+Co' 26-^^/2) _ 2g^^M,(l - 4pe-P) + I^^M.j j j | (/.(n, q^) |2, (2.3) 

whereto ^ £o+2r-HCo/c^)^ (l-2e-P/2), Hereafter, we assume the relation F > 2^q / d? , promisingH^ 
the presence of the p > 1 region in which the vortices tend to be confined between two neighboring 
layers. Since only the LLL modes were taken into account here, the gauge field 5Ay in x-y plane 
perpendicular to B, as in B || c case, need not appear as the gauge- invariant combination Qy — 
2-K6Ay/(j)Q in 7^2- However, in the present B _L c case of the LD model, the Qy\^{n,q^±)\'^ and 
qy6Ay\ip(n,c[±)\'^ terms appear within LLL even for a smaller p as a consequence of the partial 
breaking of the LLL degeneracy due to the layering. In particular, the latter term implies that 
the linear responses to 5 Ay, at least at low enough temperatures, can be defined within the LLL 
modes. Further, it is important to note that the last three terms of eq.(2.3) are summarized in the 
form £,Q{qy — 2Tr6Ay/(j)Q)'^, just as the second term of (2.3), if the 0(e~^/^) corrections in eq.(2.2) 
are neglected by assuming a strong field or a large enough anisotropy. This recovery of the isotropy 
in y-z plane in high p should be expected in the LD model: Since, in this model, the pair-field is 
defined only on the discrete superconducting layers, the vortices in p > 1 are confined between two 
neighboring layers, and consequently, the field-induced anisotropy on the layers tends to be lost 
with increasing p. This recovery of the y-z isotropy and the faclH' that the high field limit of the 
mean field transition line given by = is the same as the exact resultS-* suggest that the high 
field regime of the LD model may be well approximated by the present LLL approximation even 
in examining the response properties in y-z plane. For these reasons, for later convenience, we will 
give the high field approximation of eq.(1.2) including the term in the formi^ 



Q \ t^=y,z 



^dLy L 2. 



X E Vb(ni - n3,n2 - 71-3, gyj) 99*(ni,q_L,i) V3*(n2,q±,2) 
Qi,Q2,Q3 

X9?(n3, q^^a) ip{ni + n2 - na, q_L,i + q±,2 - q±,3), (2.4) 



where 



V j=l,2,3 
^v{Qi - Qs) v(Q2 - Qs) E ^(Qi + ^2 - 2dm/ rl) (2.5) 
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is a bare vertex function with v{Qi) = exp(— p(nj + qy^ir'^/d)'^ /2), and Qj = Qiy + Qz^iZ = 
{dni/r%)y + q±^i. 

The above analysis is not modified by replacing the constant 6A assumed above with a y- 
dependent 6Az and/or a z-dependent 5 Ay. This case is relevant to examining the tilt response 
near and above the glass transition and will be discussed in §3. Rather, we will comment here 
on the superconducting part (/x = y or z) of dc (uniform) conductivities in the liquid region 

above where the pinning disorder term, eq.(1.3), may be neglected. According to the Kubo 
formula, it takes the form (Tg^^^ = — 9T^(q = 0,i^})/dO,\fi^^Q, where T^(ir2) is a dynamical helic- 
ity modulus defined by substituting a r-dependent gauge disturbance 6A^{t) into a quantum GL 
action^) consistent with eq.(1.2), r the imaginary time, and is an external Matsubara frequency. 
If the gauge field 6 A is spatially uniform, the derivation of a LLL quantum action with r-dependent 
5 A is the same as the above derivation of eq.(2.4). Then, since the gradient and gauge-dependent 
terms when <C 1 are isotropic in y-z plane, and the q^^i dependence {I = 1, 2, or 3) of the bare 
vertex function Vq of eq.(2.4) is negligible, it is clear that the resulting (7s,yy and as^zz in such high 
fields are the same as each other, implying the in-plane "force- free" behavior mentioned in §1. 

As is clear from the above derivation of eq.(2.4), this "force- free" behavior becomes more accurate 
with increasing p = 2TTBd'^Vf/(po, i.e., with increasing B and/or T. This is consistent with the 
experimental factilli that, in the same field range, the in-plane "force-free" behavior is seen not 
in less anisotropic YBCO but only in strongly anisotropic BSCCO. In the liquid regime in p ^ 1, 
the in-plane conductivities (/U = y, z) are ^-independent and well approximated by the 

fiuctuation conductivity for 2 D and B = case, obtained by neglecting the Kosterlitz-Thouless 
critical behavior, and have the form 

= (- J) ) ""^ (l - S) ) • P-") 

where Rg = irh/le^ = 6.45(kr2) is the resistance quantum, is a constant of order unity, and a 
relation (2.21) (to be introduced later) was used here. Namely, the resistivity in the liquid regime 
in p ^ 1 yields the thermally activated behavior with no pinning effect. 

In a previous paper,© one of us has ascribed the "force free" behavior to an extremely strong 
thermal fiuctuation (equivalent to the vortex loop excitations) enhanced by the strong anisotropy 
of BSCCO, and has argued the discrete layer structure not to be essential to this phenomenon. The 
feature suggested in Fig. 2 of ref.ll that the B dependence is lost with increasing B indicates that 
the present explanation requiring the layer structure is more reasonable, and the assumption of an 
extremely strong fluctuation, which cannot be described within the LLL approximation, seems to 
be unnecessary in understanding this resistive behavior. 

We return to detailing our theory. The isotropic response in y-z plane mentioned above, i.e., the 
"force free" behavior does not change by including any point disorder unless this static disorder 
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makes the supercurrent within the superconducting planes anisotropically random. Such a unusual 
possibility of disorder will not be considered here, and we will focus on eq.(1.3) as a reasonable 
model of such an isotropic point disorder. After replicating the model (1.2) with the term (1.3), 
the replicated effective hamiltonian0) becomes HcS = Ea=i ^ld(93^''^) + Ea,b=i Wran((/'^"\ V'^''''), 
where 

_A(2) 

Wran('/^^''\V'^*^) = r^.j j ^o("'l - n3,n2 - ng, Qy^j) 

2kBi (ILyLz o rT o 

x^(")*(ni,qx,i) ^('')*(n2, qx,2) V^"^(?^3, q±,3) ^^''^(^i + - n3,qx,i + q^^s - q±,3)- (2.7) 

Using this replicated hamiltonian, the glass susceptibility within LLL consistent with Ties is written, 
after using the Poisson summation formula, in the form 

XG(k)^AA Y: e^'^^'^(— ^) < v^Q, v^^^ > < ^Q,+k, V^^j,+k, >, (2.8) 
Qi,Q2 

where k_|_ = kyi/ + fc^i, and AA is a normalization factor weakly dependent on ky. In obtaining 
eq.(2.8), we assumed again <^ 1. 

Following the procedures used in B || c casejllil' let us consider the perturbation series of eq.(2.8), 
in particular its irreducible part In the limit of weak disorder, /ir consists of a single impurity 
line, carrying a strength A, with vertex corrections determined only by the interaction (quartic) 
term of eq.(2.4). Essential features to /jr can be seen in the studyil'S^ of xg(0) in the simplest 2 
D limit of B || c || z which will be reviewed here. In this case, /jr is a constant so that xg(0) is the 
simple geometrical series, xg(0) = 1 + /ir + Z;^ + • • •. In a representation, corresponding to eq.(2.4), 
of hamiltonian in this B || c case, /ir obtained in ref.22 is expressed in the form 

XR Ae ff 



hr\2d,± = y2 Y ^kyVp-p' (dp^Q - 2xRV{ky\p)) (^,0 " "^XrV {ky\p')) , (2.9) 



where V{ky\p) = N~^J2kx Vke~'^^^^s, k is here a wave vector perpendicular to B, is a fully 
renormalized vertex part corresponding to the bare one Vk^Vky with Vk^ = exp(— /c^r^/2), Aeff = 
A/eg^(T), with Sq* (T) = Eq^T/Tcq, is the effective pinning strength,0'ii^ and xr is an interaction 
strength proportional to the squared pair- field propagators .0' According to the mean field result 
and a more general expression of the Abrikosov factor /Sa, one notice J^'il) that, in the limit of the 
perfect (square) vortex lattice, V^^ reduces to (2x_r)~^(1 — NvJ2G^=o^k,G)j i-^., 

Vikylp) - (2x«)-i(Vo + 4„o-E^fc..c?2Ee"'"''^'^)' (2-10) 

where G = (Gi,G2) is a reciprocal lattice vector, and the total number of vortices. Eq.(2.9) 
with eq.(2.10) simply becomes /ir|2d,± esc Nv{(3a — 1). Since the factor /5a — I is always positive, 
reflecting a spatially varying 1-01 in any vortex state, this prefactor in /ir indicates that, at the 
mean field level, a glass ordering occurs as soon as the vortices solidify on cooling. Namely, a 
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spatial variation of \ip\ (i.e., the fact (3a > 1) in ordered phases, due to the presence of vortices, 
induces a glass ordering in nonzero fields contrary to in zero field case (see §1). In a more pinning- 
disordered case, the factor in the above relation is replaced with the number of vortices in a 
correlation area, A'cor, to be realized below a (if any) solidification transition (or crossover) line. 
This interpretation was justifiedii) noting the fact that a term in Vj^^^ corresponding to the last 
term of eq.(2.10) is nothing but the structure factor measuring the vortex-positional ordering. In 
a dirty enough case, the above-mentioned term in Vj^^^, just like above the melting line, becomes a 
continuous function of |k^| decaying with increasing |k_|_| even below the pinning- free melting line, 
and the positional correlation length shrinks to a value of the order of (Note that, as far as the 
fluctuation is weak or negligible, a nonzero positional correlation length is well-defined). Then, the 
factor A^cor is foundii) to decrease to a value of 0(1), resulting in a decrease of the glass transition 
point. It is this situation with A'cor of order unity which is accidentally describedil^ in the simplest 
ladder approximation with no 

i) vertex corrections to the impurity line. Further, a decrease of 
temperature, i.e., an increase of B, along the melting line leads to an increase of Acfr and hence, 
corresponds to an enhancement of microscopic pinning. The B-T phase diagram thus obtained0) 
in B II c is consistent with recent data@) of twin-free YBCO (see §4). This example for B || c 
case indicates that an inclusion of a vertex correction, implying the vortex positional correlation, 
to the impurity line organizing XG is indispensable to understanding a field dependence of the glass 
transition temperature for a fixed pinning strength. 

Now, to perform the corresponding analysis for the present B _L c case, let us first examine the 
fully-renormalized vertex part in the ideal limit of the perfect solid in p > 1 which was called the 
w = 1 regular solid in Ref.5. The mean field result of Pa of this solid takes the formS 

(^a'^ = E(-l)""^o(m,n, qyj = 0). (2.11) 

m,n 

In the cases including thermal fluctuations, Pa is expressed as /(O) + 1, where /(R) is the density- 
density correlation function, 

0) 

and takes the form& in B _L c 

/(R) = N^\J2 E G(q^) E ^o(ni - na, n2 - - M, qyj = 0) 

" <ii Qi,Q2,Q3 

x( < (/?*(ni,q_L,i) '/3*(n2,q±,2) '/3(n3,q±,3) v{ni + n2 - n3,q_L,i + q±,2 - q±,3) > 

- < lv^(n3,q±,3)P X |^(n2,q±,2)P > Sq,,q,) e'^-^^^^-^^^^ , (2.12) 

where R = {Md, Y) with an integer M. Noting that, due to the partial LLL degeneracy the 
fully-renormalized interaction term takes the form 

zai^yl^z r\ r\ n> 
Qi,Q2,<43 

X(p{ni + n2 - n3,q^,i + q_L,2 - q±,3), (2.13) 
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and writingp) eq.(2.12) in terms of this renormalized vertex V{m,n, q_Lj) and comparing it with 
eq.(2.11), one finds that eq.(2.11) is recovered from eq.(2.12) when 

-^^G(q^)G(q^ +kx)F(ni,n2, q±,ql,kx) (-1)^+"^'^^ + <5„„o + <5n„0, (2.14) 

in the hmit of vanishing q^ and k^. Notice that eq.(2.14) corresponds to eq.(2.10) in B || c. By 
noting that the renormahzed version of eq.(2.7) simply becomes 

_A(2) 

2kBTdLyL, Q„Q„Q, n„n, 
xCv{n3 - n5,ni - ns, q±,i,q±,3) Cvin2 - 71.4, m - 71.3, q_L,2, q±,i + q±,2 - q±,3) 
X(^(")* (ni, qx,i) (^(^^* (^2, q±,2) V'^"^ (^3, q±,3) f^^^ {rii + na - 713, q^,i + qx,2 - q±,3), (2.15) 
where 

C^{n,m, q±,i,q±,2) = Sn,o - ^ ^(q^ + q±,i)G{q± + q_L,2) Vin, m, q±, q±,i, q±,2), (2.16) 

and applying eq.(2.14) to this, the expression of eq.(2.15) simply becomes in clean limit 

_ a(2) (ni^ns) 

KZ{^.,^p) = ^, . N, (-l)("^-"^)('^^-"^^E^o(ni-n3,n,g,,,=0) 

X (-l)"("i~"3) ipl{ni,q_^^i) ip}{n2, q±,2) fa{n3, q±,3) 'Pf^ini + ^2 - m, q_L,i + q±,2 - qx.s)- (2.17) 

The condition ni 7^ 713, indicated above, in the summations of eq.(2.17) corresponds to a fact 
that the irreducible vertex /ir of XG is proportional to e~^ (oc (3^ — 1) (see the expression of re 
in eq.(2.19)) and implies that, as in B || c, spatial variations of and of the supercurrent in a 
direction perpendicular to B, which is the y-direction parallel to the superconducting planes in 
B _L c, in the expected solid state are indispensable to a vortex glass ordering. The prefactor 
Nd = Lx/d in eq.(2.17) corresponds to the prefactor A^^ in the resulting /|2d,± ( — ^vi^A — 1) ) 
in B II c case and appeared here owing to the ideal assumption such that the vertex correction to 
the impurity line can be expressed with no pinning disorder. Due to this large factor, the glass 
transition in clean limit of B _L c case also occurs, as in B || c case, as soon as the vortices solidify 
at the disorder-free melting line Tm{B). 

Physically, the A^^^-factor of eq.(2.17) implies the number in x-direction of Bragg spots, appearing 
below the melting line, at nonzero fc^'s (As indicated in eq.(4.5) of ref.5, the "Bragg peaks" at zero 
ky already appear above the melting line, merely reflecting the original layer structure, and are 
unrelated to an ordering in the vortex state). To justify this identification, one has only to notice 
that the origin of this A'fi-factor is a trivial (or, unlimited) summation of 71,4 or 715 in eq.(2.15) 
and that, in eq.(2.12), {n2 — n^jd plays, just like n/^d and n^d in eq.(2.15), a role of a coordinate 
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perpendicular to the layers, while Qi — Qs is a wave vector in y-direction parallel to the layers 
there. Since a positional correlation length just below Tm in real systems shrinks with effectively 
increasing disorder, the N^i-iacioT arising by assuming the limit of weak pinning should be replaced, 
in more realistic cases, with a A-dependent smaller factor. Diagrammatically, this implies that, 
for such a highly disordered case, the impurity lines within the vertex correction to each impurity 
line organizing the ladder of XG are not negligible any longer. Although it is not easy to find a 
A dependence of the vertex correction within the framework of the present work, we expect by 
following the corresponding analysis in B || c case that, in such a disordered system, the first term 
of r.h.s. of eq.(2.14) will become a continuous function of |ni| and \n2\ decaying with increasing |ni| 
or \n2\ (This can be recognized through a perturbative computation valid above T^)- Then, the 
and ns summations in eq.(2.15) becomes convergent, and a prefactor in the corresponding one to 
eq.(2.17) does not depend on a system-size any longer. Hereafter, we will replace A^^ in eq.(2.17) 
with a smaller number Ac (A) by imagining a moderately disordered system. Similarly to that 
in B II c case, this A^ physically corresponds to a (dimensioless) correlation length perpendicular 
to the layers. In such a moderately disordered system, the first order solidification transition at 
Tm is weakened with increasing p and without decreasing T because the rigidity essential to the 
solidification, i.e., the shear rigidity, is exponentially small,i^ and hence, the vortex system in higher 
p is likely to be more susceptible© to the pinning disorder. Namely, besides that in B || c (see the 
preceding paragraph), this mechanism also contributes to disordering the solid due to an increase 
of B. Hence, the first order solidification line may terminate at a p- value pc before entering the 
p ^ 1 region. In a dirtier system with pc < 1, the layering effect in p < is a small correction!) 
to the solidification, and thus, phenomena near Tm in such p's will be similar to those in B || c. 
It is straightforward to derive an expression of xg{^) in terms of eq.(2.17), which takes the form 

XG(k) = /q^ Gq,Gq^+k, (1 " /ir(k) with 

/i,(k) = ^Ae / Gq,Gq^+k, E E(-ir"^^°^ K %.)e^'^"'=^. (2.18) 

Here Gq^ = ( Ai + Co ^1 )^^ is the propagator in p » 1 of the pair-field fluctuation in LLL. Keeping 
the lowest order terms in e~^ in the n, m-summations of eq.(2.18), the glass susceptibility is given 
by 

XG(k) a 1^ re + l^ki + 1 - cos ^2fc,d j j , (2.19) 

where = a — 1 with = AAcexp(— 2p)/(27r/u). If we focus on the p-region in which, due 
to point disorder, a first order transition!) accompanied by both the positional and (ordinary) 
superconducting orderings does not occur at — Tg^ (see below) any longer, the renormalized 
mass ^ must remain positive (noncritical) even in T < Tm because, as reflected in the expression 
(2.4), the LLL superconducting fluctuation is two-dimensional in nature. Owing to eq.(2.14), the 
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mass-renormalization can be performed in the same manner as in ref.22, and one easily finds 

^^ = ^^o + 2ne'i\T) {^e-P^'f ( - Aeff {Ha - 1) ) ■ (2.20) 

Since — 1 ~ 0(e~*') ^ 1 due to the vortex-confinement between neighboring layers, wc can 
neglect the last term of cq.(2.20) even at low enough temperatures so that the solution becomes 
A-independent. The resulting ji can be expressed below the mean field phase boundary in the form 

Here /io = /^o + £g^(^) 1ii(2a*c/£g''(^))/^' ^^"^ /^c is a cutoff constant of 0(1). 

Using this /x-expression, the uniform glass susceptibility XG(k = 0) is found to diverge at the 
Josephson- vortex-glass (JG) transition temperature Tjg{B), where 

^H^ = 2.sg)^^(5- 4-)(A)), (2.22) 

where Eq^ = 16Tr'^X'^{0)kBTco/{<pQd) is the 2 D fiuctuation strength defined in §1, and Tc2 is the 
mean field transition line. For simplicity, the correlation volume Nc included in the characterstic 
field 

5f ) (A) c '^V ln ( MA^] , (2.23) 

which is the only measure of a disorder strength in eq.(2.22), is assumed here not to be S-dependent 
(see, however, §4). In very clean systems, Bc^^ becomes quite large, while it decreases with in- 
creasing disorder and may become negative in moderately disordered case because the inequality, 
a < 1, is expected even in real systems with strong thermal fluctuation such as BSCCO. 
As clearly seen in the above discussion, eq.(2.22) is well defined in S > Bc^^ and in Tjq < T^. 
The JG transition temperature Tjq approaches the melting line as B decreases and approaches 
B^^\ or equivalently, as the disorder becomes weaker at a fixed B value, i.e., when grows and 
approaches N^i- Oppositely, an increase of B leads to a decrease of iVc, and hence, Tjg{B) goes 
away from T^. In high enough fields, Nc{^) reduces to a constant of order unity as in B || c, and 
eq.(2.22) is simplified in the form 



TjGiB) l-£g)lnVA72^ ^0 



where Tc2 = 7co/(l — Eq In^ A/2tt). Note the weak dependence of Tjg on the pinning strength A. 

It is valuable to point out that, except a difference in the numerical prefactor of order unity mul- 
tiplying B, eq.(2.22) is of the same form as a B-T line following from the disorder-free Lindemann 
criterion 

4 

< 4 >har~ ^, (2.25) 
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if Br' in eq.(2.21) is replaced by B^^> ~ -{(l)o/{2TTd^Vf) )ln[eV(|/iok^)]- Here, Sy = r^d^x (see 
ref.26) is the displacement fluctuation parallel to the layers, and x is the harmonic phase fluctuation 
around the w = 1 regular solid, i.e., the disorder- free ground state in p > 1, with its fluctuation 
energy 

^^(.=1) ^ k^TM /dvf 5: {d,xf + ^o'\f^o\e-^{d'd'^xf] (2.26) 

(see eq.(3.8) of ref.5). The thermal average < >har in eq.(2.25) is defined in terms of this 

fluctuation energy. Further, we used in eq.(2.25) the fact that the average spacing between the 
vortices in y direction of the w = 1 solid is 27rr^/d. 

To simplify our discussion on the result (2.22), we assume for a moment S^'"^ <C (t)o/ {2TTd'^\/f e%'') 
(i.e., e^hn\/A/27r < 1) and hence, identify eq.(2.22) with 

r,2-Tei(p) ~Tei(p)eg)p. (2.27) 

This Tci{B) decreases and approaches zero with increasing p. On the other hand, the supercon- 
ducting transition in each layer in the disorder-free case may occur at a nozero temperature even in 
large enough p (oc \/r) because the 2 D Meissner response is permittedJ^l' Thus, the disorder-free 
superconducting transition point Tsc should exist far above Tei, at least in large enough p, and the 
only possible parameter-dependence determining such a Tsc will be 

Tc2-T,,^T,,e'§\ (2.28) 

Further, the solidification (melting) transition at Tm was argued in ref.5, as in B || c case, to be 
of first order and simultaneously the superconducting transition, implying that = Isc as far as 
the thermally excited vortex loops are negligible (see Note added in poof in ref.5). Thus, if the first 
order freezing transition line at Tm terminates, for instance, due to the point disorder, somewhere 
in the B-T phase diagram of a sample, the ordinary superconducting ordering in B _L c will not 
occur at any nonzero temperature in this sample. 

However, a possibility (in the disorder-free case) of a splitting between the superconducting 
and melting transitions (i.e., < Tsc) was not excluded there, at least if going beyond the 
LLL approximation. If Tm < Tgc, this melting transition, occuring between two superconducting 
phases, is associated onlyi^ with the positional ordering across the layers, and Tm will decrease with 
increasing p. A reasonable guess of such a Tm will be that Tm oc Tei(p). However, it was arguedil^ 
that a B-T line obtained in terms of a Lindemann criterion neglecting vortex displacements across 
the layers should be far below a true melting transition in B _L c. Actually, such a possibility 
that Tm — Tci(p) < Tsc in p > 1 is unreasonable even from another point of view. To explain 
this, let us first point out that the expression (2.27) on T^i is independent of the starting model 
and valid in the phase-only approximation, as far as the inter layer shear modulus is small likeil' 
~ e~P, because the only model-dependence appears, in such a harmonic fluctuation energy, as a 
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prefactor of the interlayer shear energy, while such a prefactor is reflected merely logarithmically 
in a quantity corresponding to -B^™^ or Bc^^ and hence, is negligible for large enough p. On the 
other hand, the melting line Tm {B) in p ^ 1, where the layering effect may be negligible, will be 
well approximated by the corresponding curve of the anisotropic GL mode® in B _L c, which is, 
in the phase-only approximation, given by 

Tc2 - Tj^^ (P) ^ (P) 4'^ Vp, (2-29) 

except a numerical prefactor of order unity of r.h.s. (see eq.(2.18) of ref.26). Thus, as shown in Fig.l, 
just three curves expressed by eqs.(2.27) to (2.29) are expected in the phase-only approximation of 
the LD model. Note that the point p ~ 1, at which the three curves merge with each other, is a kind 
of dimensional (3 D to 2 D) crossover of the melting line in B _L c and corresponds to the fieldl^) 
^dc ~ in B II c case. Such a crossover to the 2 D regime in B || c is expected to occur 

with increasing T in order to avoid a unlimited enhancement of thermal fluctuation (Note that 
the B = fluctuation strength \J Eq^ oc Eq^Vt in 3 D grows unlimitedly as F increases). Besides 
this, in the present B _L c case, we have not only such a quantitative field-induced enhancemnet of 
fluctuation but also a qualitative field-induced reduction of the fluctuation due to the (1 D to 2 D) 
rise& of its dimensionality. Thus, it is quite unreasonable to expect that in B _L c and p ^ 1 
would be given by Toi(p) far below (the extrapolated) Tm^ . 

Now, it will not be difficult any longer to understand why the vortex-glass transition position 
Tjg[p) in p ^ 1 is related to the pinning-free Lindemann criterion (2.25). First, the small factor 
~ e"'' is common both to the interlayer shear modulus in the elastic energy and the prefactor in 
/ir of the glass susceptibility (see eq.(2.18)), and its origin is a weak spatial variation of {ipl in y- 
direction, i.e., a (small but) positive /?^ — 1, in the solid state. The relation /3a — 1 ~ e"^ may remain 
unchanged as far as the disorder is weak so that a granular structure, i.e., a spatial variation of 
induced by A at low enough T, on the layers is invisible. Further, the rise, due to the layering, of the 
dimensionality of superconducting fiuctuation is also reflected in the dispersion of the elastic energy, 
and the resulting exponential T-dependence, combining with the above-mentioned exponential p- 
dependence, makes a role of disorder in the resulting transition line weak logarithmically. The 
only effect of the point disorder on the vortex elasticity near and above Tjq (i.e., in a vortex- 
slush regimeBil') is regarded as a renormalization on the interlayer shear modulus oc e~'P. When 
combined with the model-independence of the Tei(p)-line, this correspondence between Tjg{p) and 

(2) 

the elastic temperature Tei(p) convincingly suggests that the behavior Tc2 — Tjg{B) ~ Tjg{B) Eq p 
may be generally valid irrespective of our use of the LLL approximation. 

§3. Responses Quantities near Tjg 

In this section, we examine linear response properties near and above the transition point Tjq- 
As mentioned in relation to eq.(2.4), the physical quantities in the vortex liquid region with e~^ <C 1 
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are insensitive to the field B, reflecting the confinement of vortices in the spacings between two 
neighboring layers in p — > oo, while in the vortex-slush regionlllil' above Tjg but below they 
should depend on B due to the remarkable S-dependence of Tjg{B) given by eq.(2.22). 

This picture on the ^-dependence also applies to the in-plane resistivities p^^ {fj, = y or z) = 
^/ {(^n + <7s,^n)i where cr^ is the normal part (quasiparticle contribution) of the in-plane conductivity 
and reasonably assumed to be ;U-independent (i.e., isotropic in y-z plane). Since, just as in the 
case with no disorder (see the discussion leading to eq.(2.6)), the gauge-coupling and the gradient 
terms in e~*' ^ 1 is approximated in an isotropic form in y-z plane, and the dependence of the 
hare four-point vertices in the replicated hamiltonian in such high fields will merely play secondary 
roles, the total conductivities in high fields B _L c are also isotropic in y-z plane. Namely, we have 
the "force free" behavior in e~^ <C 1 even when the point disorder is not negligible. Further, the 
recovery of the gauge-coupling form in y-z plane (see eq.(2.4)) within LLL also implies that this 
gauge-invariant gradient operates on the glass order parameter formed within LLL. This enables!^ 
us to, just like the conductivity!) parallel to B in B || c case, anticipate the diverging behaviors 
of near Tjq according to the scaling argument same as in the zero field transition, and 

consequently, we expect the critical behavior 

as,^^r.{T-TjG{B)r'^^-''\ (3.1) 

where zj{> 4) and vj{> 1/2) are, respectively, the dynamical exponent and the exponent of the 
correlation length, ^j(T) ~ (T — Tjg{B))~'^-' ^ describing the JG critical behavior. In p ^ 1, the 
behavior (2.6) far above Tjg (but below Tco) is expected to, on cooling, smoothly change close to 
Tjg into the behavior (3.1) which is also isotropic in y-z plane. Note, however, that, due to the 
-B-dependence of Tjg{B), the in-plane resistivity should show a B-dependent "force free" behavior 
near Tjg- 

We will focus in the remainder of this section on the static responses to a tilt of the applied 
field which are measures of the presence or absence of a transverse Meissner effect (TME) in the 
resulting glass (JG) phase. This static tilt response is a key quantity for distinguishingS^ between 
various glass phases and is in this sense more nontrivial than the resistivities which vanish at a 
glass transition irrespective of the current direction. The static tilt response AC44^^x, corresponding 
to the single vortex part of a tilt modulus, in the x-direction is definedlS by 



AC44 = b'^ '^^^^^ — ky — 0, kz 



h.2 



(3.2) 



where Ty{kx = ky = 0, kz) is the helicity modulus defined by a gauge disturbance 5Ay{z) dependent 
only on z, and the corresponding AC44^y is similarly defined in terms of 5Ax{z). Another kind of 
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tilt moduli are defined in terms of 6Az in the fornS^ 

j2 '^z i^x — ky, kz = 0) 



c: 



U2 



(3.3) 



which is a response to 6Az{y). In the disorder-free liquid regime in B || c, AC44^^ coincided^ with 
6*44^^, and AC^^^x in the (disorder-free) liquid regime in p » 1 is expected to, just as CTg^^^ (see 
eq.(2.6)), grow exponentially on cooling, i.e.. 



AC44,x ~ exp(rco(egVVr)- 



(3.4) 



First, we will anticipate a critical behavior of AC/^/^^x with no detailed calculation and simply 
by comparing with 6*44^. To do this, one only has to note that the derivation of the quadratic 
terms of eq.(2.4), in which 5^4^ (/i = y and z) are constant, is trivially extended to the case with a 
z-dependent bAy and a y-dependent bA^ ~ 6a{qy)e"^y'^ with \qy\ <^ d/r\. Namely, in this case the 
quadratic term of eq.(2.4) is replaced with 

2 

+ /^o|¥'(n,r±)|^ 



Wld,2(95) = j dyd^^^^o ^"i^y - ^'^"^^^^^^j 'fi^^^^i-) 



Hi 



-id. 



2tt 



5Az{y) (p{n,rj_] 



(3.5) 



where the real space representation for (^-fields is defined here by assuming any y-dependence to have 
wavelengths longer than r^/d. As explained elsewhere,!^ when a vortex-glass transition will occur 
within LLL, the presence of a gauge-invariant gradient —id^ — {27: / (j)Q)6A^ operating on the LLL 
pair-field in a hamiltonian such as eq.(2.4) or (3.5) implies that the corresponding gauge-invariant 
gradient operating on the glass order parameter must appear in a resulting effective hamiltonian 
TLeE, G expressed only by the glass order parameter, because the analysis in LLL leading to TLeS, G 
is formally the same as the corresponding one in zero field case. Since this similarity to the zero 
field normal-Meissner transition implies that the scaling argument© on the critical behavior of 
the conductivity cJs^^^ near the glass transition is applicable just as in zero field case, we have 
concluded^ that the helicity modulus T z{c{ = 0) parallel to B should be finite in any vortex-glass 
phase. Since, in the present case, this applies also to the y-direction due to the recovery of the 
isotropy in y-z plane, we can conclude that, as in any pinned Josephson vortex solid with no 
disorder Ty(q = 0) should be also finite in the JG phase. Based on eq.(3.2), this directly implies 
that the JG phase should have a TME perpendicular to the layers, and consistently that the tilt 
response (i.e., a diamagnetic susceptibilityil) ) AC44^a; should diverge, just as in the normal-Meissner 
transition, in proportional to the correlation length ^jg{T) on approaching Tjq from above. 



AC44,. a tjG{T) ~ (T - TjG{B)r 



(3.6) 
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which will reduce to eq.(3.4) far above Tjq{B) in p ^ 1. 

Next, let us turn to AC44_y. In this case, it is convenient to return to the original model (1.2) with 
no 8Afj^ (/X = y, z) and with replacement ij^j — > ij^j exp(i(i27rj 5Ax{z) / 4>q). After gauge-transforming 
il^ix) e'x.'piyi-Kx 5Ax{z) / (po) — > iIj{x) in the second representation of eq.(1.2) and, for simplicity, as- 
suming 6By = dz6Ax to be constant, we have 



/+00 
dV J2 e^^^'"^/'^ 




eo\m\'+eo E 

2 



2 

27r 

-id^ + —{By-6Byz)x\tP{r) 



|V'(r)-^(r + dx)|2 + ^| 



(3.7) 



in place of the second representation of eq.(1.2). When the disorder term (1.3) is included in 
addition to eq.(3.7), the inclusion of 5By merely implies a rotation of the applied field direction 
within y-z plane in which the original system in zero field is isotropic. Namely, since this disturbance 
is the same as a tilt in an isotropic 3 D system with point disorder, no TME is expected for a tilt 
within the layers, i.e., ACM,y is nondivergent at and below Tjg- Note that this conclusion was 
obtained without using the LLL approximation and hence, should be valid at any p-value. 

In contrast to this nondivergent AC44^y, the above conclusion on the divergent AC44^^x was ob- 
tained in terms of the high p approximation. Hence, one may suspect that a vortex-glass phase in 
lower p- values, where the layering is quantitatively negligible at least above Tm, will have no TME 
as well as the ordinary vortex-glass (VG) phase. 

0) 

However, as mentioned below eq.(2.2), the cross 
term ~ —Lp*[6Aydyip, which determines the current vertices in the response function, exists within 
LLL even without assuming a large p. This rather suggests the presence, irrespective of the p- value, 
of TME for a tilt perpendicular to the layers in the ordered phases. 

To clarify these points, we examine the tilt responses in the uniaxially periodic GL mode® with 
point disorder 

Wper = Wgl - J d^r[upCos{2Trx/d) - u,(r)]|V(r)p, (3.8) 

where "Hgl is the ordinary isotropic GL model, the positive constant Up is the strength of a periodic 
variation of Tco playing similar roles to the layer structure of the LD model, and Ur is a random 



potential expressing the point disorder and has the ensemble defined by = and Mr-(ri) Wr(r2) = 
^QA6^^\ri — r2). Hereafter, we assume a weaker magnetic field satisfying rs > d so that the 
ground state in a fixed field becomes a w > 1 pinned solid or the floating solid. i'il' Note that, in 
this model, the pair-field can be nonzero not only on the "layers" but also between two neighboring 
"layers" and hence that the point disorder existing between the layers is also effective as pinning 
sites enabling the vortex lines to locally deviate from directions parallel to the layers. By contrast, 
in the LD model the point disorder, as well as the pair-field, is restricted on the layers and hence, 
will not play any role for releasing the vortices locked in the layers. Due to this difference, described 



16 



in Fig. 2, in the disorder configurations between the two models, one may expect that the point 
defects existing between the neighboring "layers" in the periodic GL model compete with the layer 
structure inducing^^ the nonvanishing helicity modulus Ty(q = 0), implying TME perpendicular to 
the layers, and that this TME perpendicular to the layers may be destroyed by such defects absent 
in the LD model. Namely, a TME should occur more easily in the LD model than in the periodic 
GL model. Keeping this in mind, let us examine AC44^^ of the model (3.8) in a manner similar to 
in ref.4. 

As in previous works0& on glass transitions in B || c, we will use a high field approximation 
convenient in examining electro- magnetic responses. When treating TCql by expressing ip in terms 
of the Landau levels (LLs), the vertex on an external current perpendicular to B is accompanied 
by the next lowest LL (NLL) even if the thermodynamics is described primarily by the LLL modes. 
Since this NLL mode is not associated with any ordering in vortex systemsjll© it remains massive 
and heavy in high fields in the sense that its time scale becomes much shorter than those of LLL 
modes. In the ideal case with no source of vortex pinnings, this short time scale of NLL inevitably 
becomes the time scale for all terms of the perturbation series of a conductivity perpendicular to 
B, and consequently the vortex flow conductivity follows® This fact is clearly seen0) particularly 
in the high field approximation where the LLL modes are assumed not to interact, through the \ip\^ 
term of TCql, with higher LL modes. However, when either a periodic or a random pinning term is 
present as in eq.(3.8), a much longer time scale than that of NLL can be picked up from a section, 
composed only of the LLL modes, in Feynman diagrams contributing to a response quantity, and 
hence the presence of the time scale of NLL can be neglected in such diagrams. In refs.l and 4, 
such diagrams with a long time scale growing on cooling due to a LLL vertex correction expressing 
the VG fluctuation were found for a conductivity perpendicular to B in the case with no periodic 
pinning effect. Its examples are given in Fig. 3 (a), where the open square denotes the above- 
mentioned LLL vertex correction corresponding to the glass susceptibility. Situation is similar in 
static responses such as AC44^^, and, just like the neglect of dynamics of the NLL modes in such 
diagrams of the conductivities, the z-dependence of the NLL fluctuation propagator can also be 
disregarded in examining AC44^^. 

Once the periodic potential tx Up is included, the linear responses become anisotropic in x-y 
plane. When the gauge disturbance 5A couples to the pair-field through the GL hamiltonian, the 
tilt response AC^^^x is generally expressed in the formS'i) 

2kBT ( d\^^ 



+ < (Pl,K+k,fo,K'+k.'^hK"^o,K >]' (3-9) 

where V is the system volume. The last (second) term of eq.(3.9) vanishes in the case isotropic 
in x-y plane. For the present purpose, we need not to examine eq.(3.9) in details but only to 
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focus on how the NLL fluctuation couples to a LLL fluctuation, in terms of the potential terms in 
(3.8), outside the LLL vertex correction implying the glass susceptibility in Feynman diagrams (see 
Fig. 3). Let us consider the diagrams in Fig. 3 (b), where the NLL lines couple to LLL lines only 
through the periodic potential, by expressing the pair-field in terms of the same representation of 
LLs as in the LD model (see the sentences prior to eq.(2.1)). The NLL-LLL coupling term, marked 
with dashed line in Fig. 3, is written as 

-V2up—^ exp{-{TrrB/df) '^sm\—^Qy\ {Vi{n,<i±) f{n,q_L) + c. c. ), (3.10) 

where ipi is the corresponding one in NLL to ip in LLL. As already stated, the NLL fluctuation in 
Fig. 3 (b) can be assumed to be independent of z. Then, when ensemble-averaging over the NLL 
modes according to Fig. 3 (c), the resulting NLL propagator can be treated as a constant factor mul- 
tiplying the current vertex, and the resulting 0{6Ay) term (oc —Up6AySm{27rr'^qy/d)\ip{n,qy,z)\'^) 
can be regarded as a current vertex defined within LLL. Explicitly, if using the fact that the NLL 
propagator Gi{qy,qz = 0) =< \(pi{n,qy,qz = 0)p > in high enough fields can be replaced, ex- 
cept corrections in g^, by ?"^/(2^q) deep in the liquid regime, 0) the obtained dispersion of the LLL 
fluctuation and the current vertex coupled to 6Ay can be seen as those defined from the following 
quadratic term of an effective hamiltonian only on the LLL modes 



^per,2 = J dydz^ip*{n,r^) i Up exp{-{TTrB/df) i 1 



cos 



27rr| / .„ 27r 




+fio + Ciy-idz-'^6Aziy)] ]^in,r^). (3.11) 

Now, let us consider the diagrams in Fig. 3 (c) composed of the current vertex of Fig. 3 (b). In 
Fig. 3 (c), the left (right) diagram arises from the first (second) term of eq.(3.9). Further, note that, 
near Tjg, the gauge invariant gradient term in y direction has only to be kept up to quadratic 
order. Then, the argument given below eq.(3.3) can directly be used in the present case. The only 
difference from the LD case is the absence, at the quantitative level, of isotropy in y-z plane due 
to the difference in coefficients of the gradient terms. Since the strength Up is independent of other 
material parameters, AC44^2,. in the model (3.7) is also predicted to diverge in proportional to a 
correlation length near the resulting glass transition. Contrary to this, this divergence is absent in 
AC44^y, because this quantity is given by replacing the plus sign prior to the last term in eq.(3.9) by 
a minus sign, and hence, the contributions corresponding to the two terms of Fig. 3 (c) cancel with 
each other. We note that the diagrams, such as Fig. 3 (d), including just one Up-vertex coupling 
between the LLL and NLL modes do not give nonvanishing contributions because of the ruleU^ 
illustrated in Fig. 3 (e). Hence, only the sum of diagrams such as Fig. 3 (a) give a nonvanishing 
contributions related to the glass fluctuation which is easily verified to be nondivergent at the 
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transition for the same reason as in the VG caseEl In passing, we note that what is essential 
to obtaining the above results is the presence of a nonzero gauge coupling ~ i6Ay{z)(p*dyip after 
averaging over the NLL modes but is not the use of the approximation on Gi in getting eq.(3.11). 

Thus, we have shown that the divergent ACm^x, as well as the nondivergent AC44,y, is also 
realized in this periodic GL model. We emphasize that this conclusion is independent of the 
strength of the point disorder and of the strength B of the magnetic field which tunes the structure 
of the disorder-free vortex solid in B _L c through a commesurability condition!' with the layer 
structure. Namely, the glass phase in the model (3.7) should have the anisotropic TME irrespective 
of B. Following the consideration based on Fig. 2, this in turn convincingly suggests that our earlier 
prediction on the LD model that the resulting JG phase has TME only for tilts perpendicular to the 
layers is also likely to be valid for any B. Namely, the ground state will become the JG phase with 
a TME even in p (< l)-values where the disorder- free ground state may be a unpinned (floating) 
vortex solid, as far as the positional order in the resulting JG phase at such p- values is short-ranged. 
Further, we note that the presence of the TME perpendicular to the layers below Tjg for any p 
implies the presence in any B of the so-called lock-in transition due to a tilting of the applied field 
in real systems with point disorder. Note that the so-called lock-in transition close to B _L c was 
discussed so far just in the pinning-free case. Since the floating solid possible in the disorder-free 
case has no TME, the above result implies that the point disorder destroying the positional order 
of this solid simultaneously creates the anisotropic TME in the real JG phase. 

§4. B _L c Phase diagram and Related Issues 

The resulting B _L c phase diagrams in this paper are sketched in Fig. 4. The main point in 
our results is that the resulting Tjg must decrease with increasing B even in high fields satisfying 
B > (/)o/\/rd^, while the first order melting temperature and phenomena in the vortex liquid 
regime above Tm{B) in such high fields will be almost independenti^ of p oc -B\/f, reflecting the 
nearly complete confinement of vortices within the inter layer spacings. The origin of the field 
dependence of Tjq is a spatial variation of IV'I on the superconducting layers, which is small but 
nonvanishing in B > / {^/V d^) , and, as already mentioned, is indispensable to a glass ordering. 
Since this spatial variation strongly depends on p and disappears in high B (or high F) limit of 
the LD model, Tjg at which the resistivities vanish must depend on both B and F and approach 
zero in this limit. Namely, in B _L c case, Tjg rapidly starts to deviate from the (disorder-free) 
melting line Tm{B) with increasing B even in clean systems, and hence, the resistive curves just 
above Tjg but below also become -B-dependent. This is quite different from the (ordinary) VG 
transition lineS'© Tvg{B) in B || c case, in which the VG line Tvg{B) in higher fields (described 
as the lower dashed curves in Fig. 7 of ref.l) yields the LLL scalingEl-* similar to that of Trn{B), 
and for this reason, this TyG{B) line is often regarded, for a practical purpose, as an extrapolation 
of the Tm{B) line by neglecting a narrow windowi-* of the slush regime observed© recently in 
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YBCO. Therefore, particularly in materials with stronger fluctuation, there may be a much wider 
vortex slushSS regime in B _L c than in B || c. Further, as shown in §3, the "force free" behavior 
itself in such high fields is never affected by the point disorder isotropic on the superconducting 
layers. Consequently, the in-plane resistivity curves in p > 1 will show, near Tm, a crossover from 
a i?-independent "force free" behavior to a i?-dependent "force free" behavior below Tm- The 
B- dependent "force- free" resistive behavior below T^, is qualitatively consistent with the behavior 
at lower temperatures seen in the data of ref.ll. 

For convenience, we have distinguished a dirtier case Fig. 4 (a) from a cleaner case Fig. 4 (b) based 
on whether a critical point of the first order transition line, denoted as p = pc, satisfies pc < 1 or 
Pc > 1. First, let us consider Fig. 4(a). In this case, with increasing p, Nc may already become 
of order unity before entering the p > 1 region (Note that the shortest value of the positional 
correlation length perpendicular to the layers is the layer spacing d and hence that one can say 
iVclmin ~ 0(1)). Then, the relation B^^ < 0o/(2^(i2 ^F eg^) will be satisfied over a wide range in 
p > 1, and we will have the relation Tc2—Tjq oc Tjg B (see eq.(2.24)). Next, in turn, let us assume a 
cleaner case with pc > 1 such as Fig. 4 (b). In this case, Nc 'vap < pc will be of order N^. As already 
mentioned, this means that the system below pc is in clean limit and that Tjg — Tm there (see 
Fig.4(b)). Next, when p > pc, Nc will rapidly decrease with increasing p, because the correlation 
volume will decrease with decreasing the shear mo dulusi'© ~ e ^ in the w = 1 solid-like state. If 
naively assuming Nc ~ e"'^^^ with a constant ci > 0, the neglect in eq.(2.22) of a i?-dependence 
of Bc^^ is not valid any longer. Only in p values larger than pd, Nc will start saturating into a 
constant of order unity, and hence Bc^^ approaches a small positive or a negative constant. Namely, 
eq.(2.24) becomes valid only when p > pd- Consequently, we have conjectured the phase transition 
curves, the solid curves in Fig. 4 (a) and (b). It should be emphasized that the slope of the behavior 
BjG ~ Tco — T to be seen below Tm is gentler as the fluctuation is stronger and is insensitive to the 
pinning strength. This will become useful for examining Tjg{B) experimentally and implies that 
the vortex slush regionB'il) to be seen experimentally becomes wider in materials with stronger 
fluctuation. According to ref.l, such an extension of the temperature window of the vortex slush 
region due to an enhanced fluctuation should also occur in B || c case,i) contrary to the Bragg glass 
scenario 

i) 

of the B II c phase diagram, and, in fact, was recently observed0 through a doping 
dependence of the B-T phase diagram of YBCO. 

Theoretically, studies of a glass phase in B _L c were performed in previous workJl'l^ by assuming 
a dislocation-free glass phase. However, it was arguedS that the resulting transition is not accom- 
panied by a divergence of a glass correlation length, just like a vortex liquid-Bragg glass transition 
argued as a scenario of the flrst order transition in B || c case. Further, since the assumption in 
refs.2 and 3 that the vortices cannot move across the layers is valid only inS p > 1, the situation 
assumed in these workJ'i is, at most, limited to the narrow region 1 < p < pc of Fig.4 (b). 
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As shown in ref.5, the positional correlation of vortices perpendicular to B in B _L c first develops, 
on cooling, in y direction parallel to the layers. This is a consequence of the recovery of isotropy in 
y-z plane in high enough fields, emphasized in §2, which is essential to understanding the force-free 
behavioilllllll' of the in-plane resistivities. On the other hand, the above-mentioned feature that the 
positional ordering first grows not across but rather along the layers contradicts the conjecture!^ 
favoring a smectic liquid phase. Namely, the smectic liquid picturJ^ is incompatible with the force 
free resistive behavior. 

In ref.38, Hu and Tachiki still argue that the disorder-free superconducting transition in B _L c is 
continuous for cases with larger anisotropy in which the anisotropy-induced vortex loop excitations 
seem to play an important role in the transition, while it is reasonably of first order for smaller 
anisotropics (see Note added in proof of ref.5). It is remarkable that the transition temperature 
they found is insensitive, at a fixed field, to the anisotropy even in 1.5 < p < 2 studied in ref.38 
(see Fig. 3 there-in). This is consistent with our prediction!) that, in p > 1, Tm{p) approaches 
Tsc independent of p oc B\^, because, as commonly seen in various workJl3ilii) on the phase- 
only model, the magnetic field B always appears only as the parameter p. Hence, the transition 
temperature in ref.38 is identified with — Tsc- Clearly, this transition temperature has nothing 
to do with the observedlllii) superconducting transition temperature in BSCCO which significantly 
decreases with increasing B and also dependJll* much on the doping level, i.e., on the anisotropy 
(see below). 

After preparing the first draft of the present paper, we were aware of recent data of a.c. sus- 
ce ptibility© and the resistivities, both of which can be direct probes of a glass transition, 
in field configurations including B _L c case of clean crystals of high Tc superconductors. We will 
not discuss here complicated behavior^i* in fields tilted from the layers but focus on the data in 
B _L c relevant to the present work. In ref.l9, the superconducting transition field in B _L c was 
determined in BSCCO crystals as the position at which the lock-in behavior (i.e., the TME in 
B _L c studied in §3) seen in a.c. susceptibility data disappears with increasing the in-plane field. 

The resulting line H^'^(T) roughly obeys TH^'^ ~ (Tco — T)/T, which is the same one as eq.(2.24) 
or (2.27) if the prefactor F is replaced by \/T However, this difference in the prefactor does not 

(2) 

preclude this identification between H^"^ and Bjg because the fluctuation strength Sq proportional 
to the squared penetration depth also increases, as well as the anisotropy F, with decreasing the 

(2) (2) 

doping level. Using a resonable Eq value in under-doped BSCCO, Eq ~ 0.2, corresponding to 
the in-plane penetration depth A(T = 0) ~ 3000(A), the H^'^{T)-lme in the underdoped case 
of ref.l9 is quantitatively consistent with eq.(2.27), i.e., Bjg{T) obtained by assuming Bt^ « 
(/)o/(27rd^\/r £(7^) — 3 (T). On the other hand, the corresponding transition fields in BSCCO 
were determined in ref.39 as the positions at which the in-plane resistance at He = vanishes. 
According to ref.39, this transition is continuous. The resulting transition field in an optimally- 
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doped sample seems to be comparable with the data in ref.l9 and is roughly linear in temperature 
just like eq.(2.27). Further, in resistive measurements in YBCO,© resistivities in all directions 
were found to vanish at the same temperature Tj. This is an evidence of a continuous JG transition 
occuring below the disorder-free melting line (see §1). It seems to us that, in ref.40, the out- 
of-plane resistivity (/Jc) vanishes more rapidly (with a smaller exponent) compared with the in- 
plane resistivities. Explaining this behavior seems to be theoretically interesting and will be tried 
elsewhere. 
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Fig. 1. Three possible characteristic curves, Tf.i{p), T^'{p), and T^c, in the p v.s. T diagram of the phase-only 
approximation of the LD model. These curves merge roughly when p ~ 1. The solid curve is believed to be the 
true transition curve in the disorder-free case. 
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Fig. 2. Configurations of point disorder in (a) LD and (b) the periodic GL models. Figures are described in x-y 
plane. Solid lines, in x direction, denote the "superconducting layers" , and the open circles denote the point defects 
for the pair-field. 
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(d) (e) 



Fig. 3. Feynman diagrams on tilt moduli in the high field approximation. The dashed hne denotes the impurity 
line occuring after averaging over the point defect configurations and carries the factor A, and the connection marked 
by an open square between a NLL propagator Gi (q) (the double-soHd hne with arrow) and a LLL propagator (the 
sohd line with arrow) implies the periodic potential with strength Up. The hatched square and open sphere in (e) 
imply four-point vertices composed only of the LLL modes. See the text for further details. 
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Fig. 4. Schematic p v.s. T phase diagrams of real layered superconductors in B || layers for (a) dirtier and (b) cleaner 
cases resulting from the present study. The value pc, the p value at which the first order transition (expressed by 
the solid curve on Tm{p)) ends, measures the disorder strength. The relation (2.24) will be satisfied entirely in 
p > Pc in the dirtier case (a), while it is valid only in p > pci(> Pc) in case (b). In each figure, the left solid curve 
Tjg{p) is the JG transition line, while the right dashed curve Tc2{p) indicates the Hc2{T) curve. 
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